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Abstract 

This paper presents a model which can be used to predict the two- 
dimensional nonlinear behavior of bridged cracks in orthotropic strips. 
The results obtained using a angular integral equation formulation 
which incorporates the anisotropy rigorously show that, although the 
effects of anisotropy are significant, the nondimenaonal quantities em- 
ployed by Cox and Marshall (1991) can generate nearly universal re- 
sults (R- curves, for example) for different levels of relative anisotropy. 
The role of composite constituent properties in the behavior of bridged 
cracks is clarified in this paper. 

1. Introduction 


The bridging of cracks by fibers is an important toughening mechanism in 
fiber reinforced brittle matrix composites (BMC). Similar characteristics axe 
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shared by ductile matrix composites (DMC) with weak interfaces upon cyclic 
loading (Evans, 1991). This fact has led to a significant amount of research 
in the area of stress analysis of bridged cracks. The results calculated with 
the m odels presented by Marshall, Cox and Evans (1985), Marshall and Cox 
(1987), Cox and Marshall (1991), Cox (1992) and Ballaxini and Muju (1991) 
have provided a good understanding of the fracture mechanics of bridged 
cracks in finite specimens. In all these models , excluding the orthotropic 
finite element model used by Ballaxini and Muju, the anisotropy which is 
inherent in fiber reinforced composites was incorporated in an approximate 
way which may or may not provide an accurate account of its effects. 

The integral equation model for a bridged crack in a finite specimen de- 
veloped by Cox and coworkers can be summarized as follows. A fictitious 
line load P per unit width of crack front was placed on one crack sur- 
face at position x measured normal to the crack front. Using Castigliano’s 
theorem the crack opening displacement (COD) u(x) was determined as 
u(x) = limp , 0 ~, where W is the elastic strain energy of the system per 
unit width of crack front. This energy was then written in terms of the 
crack leng th a, strain energy release rate Q and stress intensity factor IT, 
i.e. W = JZGda' = J£K 2 /E’da'. In this equation E' is given, for an or- 
thotropic material with a plane stress crack parallel to the principal material 
axis on-axis, as 1/E 1 = yJhJ^i^/hi + (2 In + J 6 e)/21u} 1/2 , where the 
kj axe complian ce coefficients which enter the inverse Hooke’s law, as will be 
explained in the next section. Cox and Marshall made the following approxi- 
mation. They expressed the stress intensity factor K in terms of the Green s 
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function G(x, a',w) for the isotropic material, where w is a relevant char- 
acteristic leng th which introduces finite geometry. The representation K — 

2 G{x, a', w)a{x) dx, where cr(x) is the traction on the crack surface, leads 
to the integral equation u(x) = (4/J5') g da!{$ G{x\ a ', w)c{x) dx , }G{x, a 1 , w), 
which they solved numerically. This integral equation can be more directly 
derived using the weight function technique (Bueckner, 1970). Thus, the 
Green’s function G(x, a', w ) can also be referred to as the Bueckner weight 
function. This model, by neglecting the dependency of the Green’s function 
on the compliance coefficients, leads to the result that the COD is inversely 
proportional to E' . Moreover, normalization of this integral equation leads to 
the so-called bridging length scale a n , that is associated with a fully bridged 
crack in an infinitely extended material. A crack of length a << a n is con- 
sidered short in the sense that its bridging zone is still developing. For these 
short fully bridged cracks the forces in the fibers are relatively small, since 
the crack opening displacements are small. Hence the fibers do not shield 
the crack tip significantly, and the stress required to propagate the matrix 
crack is inversely proportional to the square root of the crack length, as for 
the monolithic matrix. For a >> a n the crack is termed long, because the 
bridging zone has fully developed and a steady state is reached for matrix 
cracking (the stress required to propagate the matrix crack is independent of 
crack length). 

The authors have recently learned of the significant contribution to the un- 
derstanding of the role of material orthotropy in fracture specimens for com- 
posites made by Suo et ai. (1991). They found a spatial rescaling that could 
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be used to reduce boundary value problems (BVP) involving ortbotropic 
materials to equivalent problems in materials with cubic symmetry. Under 
certain conditions the cubic-symmetric materials may be approximated by 
isotropic materials. Therefore the rescaling technique can be used in con- 
junction with existing isotropic material solutions to construct approximate 
orthotropic material solutions. These spatial rescaling relations, which will 
be discussed subsequently, were used by Bao et al. (1992) in tbeir analysis of 
commonly used fracture specimens to investigate the interplay between ma- 
terial anisotropy and finite geometry. In their finite element method analysis 
the role of orthotropy in fiber bridging phenomena was not considered. 

The work presented in this paper was initiated to investigate the effects of 
relative anis otropy on the behavior of bridged cracks in long strips so that the 
assumption used in the aforementioned analyses of bridged cracks could be as- 
sessed. Moreover, an efficient and highly accurate analytical model based on 
theory was sought that will facilitate the analysis of experiments 
conducted on beam specimens. The results obtained using a singular integral 
equation formulation which incorporates rigorously the effects of orthotropy 
show that the weight function G(x', a', w) depends not only on E but also 
other stiffness coefficients. However, the numerical calculations presented 
in this paper indicate that for the unbridged cracks in different orthotropic 
strips the CODs are indeed approximately inversely proportional to their 
values of E. A wide range of computations for the nonlinear behavior of 
bridged cracks also supports to adopt simple approximate weight functions 
which are derived for the corresponding isotropic specimen. These results 
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generalize the conclusions made by Bao et al. to fiber-bridged cracks in long 
strips. Evidently, this conclusion will facilitate the analysis of R-curves for 
bridged cracks in orthotropic beam specimens. 

2. Formulation 


Consider an infinite orthotropic plate whose principal material axes coin- 
cide with the coordinate axes. A unit dislocation with Burgers vector b x — 1 
is embedded at P(0,yo)- The stresses in the plate are (Milne-Thomson, I960): 


A 




= -| Re t^ 3 - + (i) 

where 

zi = x + Ai(y - yo), *2 = x + A 2 (y — Vo)> (2) 

Ai and A 2 axe two roots of the characteristic equation 


JiiA 4 + (2iia + lee) A 2 + In — 0 (3) 


with Im(Ai) > Im(A 2 ) > 0, A 3 = -Ai, A 4 - -A 2 . The hj axe the anisotropic 
compliance coefficients which enter the inverse Hooke’s law: 
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In equation (1), the constants Ai and A2 are determined by 

* • _ 

PlA\ -\~PzA2 — QiA.% — Q2A2 — 2/7ri, 

| 6\A\ + 62A2 — J1A1 — 72-4.2 — 

where 

Pl =ll 2 + IllAf + *(?22 + Il2Ai)/Ai, 

P2 = I12 + ^11^2 + *(?22 + ^12Al)/A2, 

= Il3 + lllAf + *(l22 + IxsAD/Ai, 

§2 = li2 + biA| + *(^22 + IiaAD/Aj, 

Si = (l+ *Ai)/2, S 2 = (l+ *A 2 )/2, 

71 = (1 — «Ax)/2, 72 = (1 ~ *A 2 )/2. 


( 5 ) 


( 6 ) 


The first equation in (5) represents the necessary displacement jump condi- 
tion, while the second represents zero net force on the dislocation. 

. When the dislocation is located in an orthotropic strip (Figure 1), an 
solution must be superposed to satisfy the traction free boundary 
conditions on the surfaces y m 0 and y = h. The stresses from the additional 
solution, denoted by superscript (2), must satisfy 

<rg>( a:,0) = -er$(x,0), 

0$ {x, h ) = («, h), cr@j (x, h ) = (x, h). (7) 


This part of the solution can be determined using Fourier transformation 
te chniques . Its displacement components are expressed as 

tt (a )(s,y) = ^ jT U&v)sm£xdx, 
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vW(x,y) =§r^-) cos £x dx , 


= EB,(e) e xp(s^ff), 

,=i 


V(& y) = ^diB^Oexp^y), 
i=i 


where 


di = [ft»? - (Aft - A J )«,]/ft, (10) 

and the Sj (j = 1,2,. . 4) are the roots of the equation 

8< + ^-ftft-l s a + A , t (11) 

P% Pi 

In (10) and (11) 

Pi = &11/G12, Pi — b&JGm Pi — 1 + &12/G12, ( 12 ) 

where the &j ; - and G12 are the stiffness coefficients of the material: 

( \ / u. K. n \ / _ \ 

CT XX Oil v cjpjg 

<Typ = 6x2 &22 0 % • (^) 

^ O'aey y \ 0 0 G12 2€*y ^ 

Using condition (7), the coefficients Bj(£) (7 = 1,2,..., 4) can be determined 
by solving the set of linear equations given in the Appendix. 

Along the y-axis the resultant stress a xm in the strip is 


cr,m{ Q,y) = + G(y,yo) 

yo-y 
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with 


G{y,yo ) 






k = —Re (A 1 A 1 + A-iXi) /%. 


(15) 


The stiffness parameter k plays a crucial role in defining the level of anisotropy 
of a given material system. After comparing it with the relevant quantities 
E' and A which are adopted by Cox and Marshall (1991) and Budiansky and 
Amazigo (1989), respectively, it is found that 


B = AEu 

4x 4x(l - i/* ) 


( 16 ) 


where u m is the matrix Poisson’s ratio, Eu the longitudinal composite Young’s 
modulus. Because the crack is now perpendicular to the x-axis, in equation 


(16) 


_1_ _ ^nl22 ^i/2^1n^/ 2 + 21i2 + h g p/ a 
E' 2 l%2 2/22 

This q uan tity was first introduced by Sih (1968). 


(17) 


With (17) the relative in-plane orthotropy is specified through three ratios: 
I22/I11, Ivt/hi and Suo (1990 a,b) introduced two parameters 


x hi _ 2ii2 + i«a 

fa’ 

as the only two parameters needed to quantify the level of orthotropy. A = 
p = 1 for isotropic solids and A = 1 for solids with cubic symmetry. Using 
Lekhnitskii’s formalism (1981) he showed that in simply-connected sheets 
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with traction prescribed on the boundary the governing equation for the 
Airy stress function can be written as 


dHJ 

dat 




( 19 ) 


This suggests that the stresses depend on material properties only through p 
and A. Furthermore, by rescaling the x-axis by f = A 1 / 4 ® the A dependence 


ran be extracted explicitly, so that in the transformed plane the stresses 
depend only on p, i.e., 


&U 
d ? 


+ 2 p 


&V 

d?dy* 



( 20 ) 


Obviously, solution to BVP for the class of materials obeying p = 1 (A 
1) can be constructed from the solution to the corresponding BVP for the 
isotropic material. It can be shown that 

<■» 

for a crack in the y-direction. The parameters A and p will be considered in 


subsequent calculations. 

Consider a cracked strip shown in Figure 1 and introduce parameters i 
and r such that y = b + at/2 and y 0 = b+ar/ 2. The discrete dislocation is 
replaced with a distribution of dislocations 

BJt) s -§£k»V) _ u (0-, t )]. (22) 


This representation enables one to write the following singular integral equa- 
tion for the traction boundary condition along the crack, with expression 
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(14) being the appropriate kernel: 

J_ t dT + f_ x K (*’ t ) b *( t ) dr = -«(*), (23) 

where K(t, r) is the regular part of the kernel 

K(t, t) = ±G(b + at/2, b + ar/2), (24) 

and the loading term is given by 

g(t) = i[o A (t)-p(t)]. (25) 

#v 

Here cta( 0 is the stress caused by the applied load in the crack-free strip and 
p(t) is the dosing stress of the bridging fibers. In general, p(t) = p{u(t)) 
where u(t) is the crack opening displacement at t. For internal cracks, the 
crack closure condition 

j^B z (r)dr = 0 (26) 

should be supplemented to equation (23). The square root stress sing ular ity 
at crack tips is modeled by expressing B x (r) as 

B,{t) = <j>(r)fVl - r 2 (27) 

for internal cracks and 

B t (r) = <f>(r)/y/l — t (28) 

for edge cracks, where 0(r) is a regular function. 

It is important to note that unlike the isotropic case, where the kernel 
of the singular integral equation is independent of the elastic constants, the 
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kernel which appears in (23) depends on the ratios EnfE\\ y Gn/En and vn 
for the plane stress case and E^sfEii, E&/Eu, Gn/En, 1/12 > P 13 and I'm for 
the plane deformation case. Nevertheless, from a viewpoint of the orthotropy 
rescaling, in the transformed fy-plane the stress field for a unit dislocation 
b m = 1 embedded in a long strip is determined by governing equation (20). 
The corresponding boundary conditions are 



&U 


= 0 


(29) 


at y — 0 and y — h. On the dislocation, the zero net force condition and the 
displacement jump condition are (Bao et al., 1992) 

[fl-0, [fk = 0, (30) 

and 

(3i) 

respectively, where the symbol [ ] L denotes the increment received on passing 
once round a closed curve enclosing the dislocation and x ^e auxiliary 


function defined by 

a 2 * 


&u &u 

d? + dy* ’ 




d£dy ~ d? * dy 2 1 ''' r '8&p 0? dy* 

Therefore, the kernel of equation (23), which is equal to (1/ K){d 2 UJdy 2 ) y 
depends on elastic constants only through the non d i me n s i o nal parameters X 
and p. Moreover, because in Figure 1 there is no characteristic length in the 
2 -direction the kernel actually is a function of p only. 
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The above integral equations are solved numerically. By approximating 
as piecewise quadratic polynomials, equations (23) and (26) are reduced 

to: 

= ~M‘<) - pM*)) (33) 

i=l 

where t, is the integral point, U is the collocation point, and the matrix 
element My consists of the weights given in Gerasoulis (1982). For edge 
cracks i « 1,2,. . .,n, and for internal cracks * = 1,2,. . .,n- 1 with the nth 

equation 

£M„ j 0(r i ) = O. (34) 

jf=l 

coming from (26). In the following discussion, we only consider edge cracked 
and center cracked strips. 

Equation (33) represents a set of nonlinear equations for the unknown 
distribution. A more efficient iteration procedure is obtained by 
deriving a complianc e matrix, as outlined, in Ballarini (1986). 

From (33) and (34), 

#Tj) = -i E «'i*M’>) - lK«fa))] 0 (35) 

K k-\ 

where N = M~ X T with Af" 1 denoting the inverse of matrix M, and T is an 
n x n matrix which transforms the coordinates from the integral points to 
the collocation points. Following Ballarini (1986), by substituting (35) into 
(28) integrating J3 *(t) from tj to 1, equation (35) is cast in the form: 

u(r ; ) = — 52 Cfrfafa) — p( u ( r *))l (J — 1, 2, . . . , n — 1) (36) 

K k=l 
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where the matrix C is known as the compliance matrix of the crack which re- 
lates the stress on the crack face with the crack opening displacement. Thus, 
a sys tem of nonlinear algebraic equations for crack opening displacement 
Ufa) is obtained. 

The total stress intensity factor which includes the fiber shielding effects 
is 

Kt = ~rfKXy/ira4>(l) 

= ^TCy/va JV„*[cr A (7*) - p{u{Tk))], (37) 

k=l 

while the so-called shielding stress intensity factor is given by 

Kf = X) (38) 

k=l 

Here q = 1 for edge cracks and rj — 1/^2 for center cracks, respectively. It 
is obvious that Kt~ K a — Kf where K A is the stress intensity factor due 
to applied load only. 

3. Results and Discussion 

Zero fiber bridging 

To check the numerical procedure calculations were first performed for 
the case of zero fiber bridging stress. The computed stress intensity factors 
a.nH crack mouth opening displacement (CMOD) for uniform tension and 
bending, for various values of relative orthotropy and crack length, are prac- 
tically the rnwift as those reported by Delale and Erdogan (1977) and Kaya 
and Erdogan (1980), and are not presented here. Selected results for plane 
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stress edge cracks under uniform tension are presented in Table 1 and 2 to 
hi g hli ght the significance of material anisotropy. 

Table 1 and 2 show the normalized stress intensity factor K/<r^yj7ca/2 
and normalized crack mouth opening displacement [u]l?/2oooh, respectively, 
as functions of the orthotropy rescaling factors A and p. It is observed that 
for a given crack length, both the quantities are independent of A and are 
weak fun ctions of p. This is because the kernel of equation (23) is a function 
of p only. As aresult, the parameter p will slightly change the proportionality 
of [ti] to \/Ej . Such effects also exist on the relation between bridging force 
and relative anisotropy. 

The normalized COD profile for a plane stress edge crack of length equal 
to one half the depth of the beam is plotted in Figure 2 as functions of 
relative orthotropy which is obtained by keeping the values of E& = -Ess, 
G 12 — £^ 2 / 2 ( 1 + 1 / 12 ), »^2 as i/ia = z /23 = 0.3 unchanged while varying the value 
of Eu only. For the cases EufEa = 1, 2, 5 and 10, EfEn = 1, 1-47, 2.43 
and 3.55, A = 1, 0.5, 0.2 and 0.1, and p = 1, 1.63, 2.77 and 4.02, respectively. 
It is seen that, though the material anisotropy has significant influence on 
the COD, the relation between [u]E , /2a 00 b and y/a is almost independent of 
the level of relative orthotropy. The largest discrepancy occurs at the crack 
mouth, but even for Eu/Eka = 10 it is less than 8% compared with the 
isotropic case. 

Initially fully bridged cracks 

Results are calculated next for bridging fibers whose strength satisfies a 
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two parameter Weibull distribution and whose sliding is resisted by a constant 
frictional stress r. The dosing stress can be approximated as (Cox and 
Marshall, 1991) 

p(u) = f£U 1 ' 2 exp(-Ut m+1 '> /2 ) y (39) 

where / is the volume fraction 
the normalized COD: 

U = «/u», 


of fibers, m is the Weibull modulus and U is 


.. S 2 fl(l - f)E„ 
“ irEfEn 


s = {*m 


rn+lv 
m + 2 ’ 


(40) 


with (a) the average fiber strength, T the gamma function, R the fiber radius, 
and 2? m , Ej and the moduli of the matrix, fiber and composite in the 
fiber direction. 

The effects of anisotropy are demonstrated through the specific example 
of two plane strain orthotropic strips: strip 1, nearly isotropic with bulk 
properties Eh = 279.6GPa, E& — Ess — 253.9GPa, Gu — 97.66GPa and 
vu as i/ 13 = = 0.3, and strip 2, with all properties equal to those of strip 


1, except as 1398GPa. Their rescaling parameters are A s 0.9163 and 
p a* 1.086 for strip 1 and A = 0.1963 and p — 3.048 for strip 2. 

For each strip three different normalized crack lengths were considered: 
Aa=3, 5, 7 for uniform tension loading where 


A a= a/a». 


TcERjl-~f)E m E 

\6rfEfEu 


(41) 
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The dimensionless strip width is H = h/a n — 10. 

In each analysis the crack length is held fixed as the loading is increased. 
By controlling the length of the bridged crack results are obtained which 
highlight the transition from stable to unstable behavior of the crack. The 
resulting physical parameters are presented in Figures 3-8. These include 
the square root of the normalized crack mouth opening displacement y/Uo = 
yju(— l)/u», square roots of the normalized total and shielding stress intensity 
factors tJKt/Kq and ^Kf/Kq as functions of normalized loading parameter 
y/Sfc, where 

K* = 2/B^2a tt /3ir, (42) 

S a = <7oo//S for uniform tension under stress a «. 

It should be noted that the noudimensional variables used in equations 
(39)-(42) were defined by Cox and Marshall (1991). They are adopted in 
the present analysis to drive home the point whether this set of nondimen- 
flirmal parameters are sufficient to predict the behavior of bridged cracks in 
composite materials which possess different levels of anisotropy. 

Notice that for a bridged Griffith crack (6 = hf 2, a/h — ► 0), the effects 
of free surfaces disappear, and the Fourier transform kernel K(t ,t) vanishes. 
Consequently, the compliance matrix C in equation (36) is independent of 
elas tic m oduli . As illustrated by Figures 3-5, where the Weibull modulus 
m = 1, the relations between Uo and S a , Kt/Ko and S a and Kf/Kq and S a 
are universal for all orthotropic strips with common values of A. 

The stability of the bridged cracks can be viewed either through the y/Uo 
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versus ,/SZ, JKtJKo versus y/Sl or yjKp/Ko versus plots. Figure 
3 shows that for relatively short cracks an increase in load is needed to 
increase the crack mouth opening displacement (CMOD). For long cracks, 
on the other hand, as the fibers start pulling out the CMOD increases even 
if the load is decreased. In other words, the instability is indicated by a 
discontinuity in the CMOD. 

Consider next the stress intensity factors. Figure 5 shows that for the 
stable short cracks, the shielding produced by the fibers varies stably, while 
for long cracks there is a sharp reduction in shielding as the fibers pull out. 
This sharp reduction corresponds to the discontinuous increase in total stress 
intensity factor in Figure 4. 

Of interest is how the normalized quantities A and H work when the 
orthotropic strips have finite geometry configurations. Figures 6-8 show the 
curves y/Uo versus -y/^, Kp/Ko versus and Kp/Ko versus yfSl for 
bridged edge cracks with m = 0 under uniform tension. From these figures it 
is seen that even in finite geometries the above normalized quantities can still 
yield nearly universal results, regardless of the level of anisotropy. Although 
they axe not presented here, the results for bending loading and for center 
cracked strips showed similar trends. 

The above results suggest that the regular kernel K(t, r) in equation (23) 
can be expressed as r)[l +Y(p,a/h, b/h)] where r) is the regular 

kernel of the corresponding isotropic cracked strip and Y (p, a/h, b/h) is a 
hi gher order perturbation term which is a function of p, a/h and b/h with 
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Y(l,a/h,b/h) = 0 and Y(p, 0,0.5) = -1. To obtain the detailed form of 
Y (p, a/h, b/h) further numerical investigation is needed and is not discussed 
in this paper. 

Here it should be emphasized that although in all the above calculations 
only the bulk moduli enter (4) and (13), for real fiber-reinforced composites 
the bulk moduli axe related to the constituent properties. There are many 
procedures published in the literature to perform this homogenization process 
(Mura, 1987). 

Prom the foregoing analysis it is seen that the length parameter o» is a key 
parame ter to define all relevant normalized quantities such as the normalized 
crack length A, normalized strip width if, normalized stress intensity factors 
Kt/Ko and K F /K 0 . Then, what role is played by the constituent properties 
in determining the value of «„? To this end, the parameter o„ is written as 

v Efi 

a " s= 16 “ a 

where 

(1 - f)E m E 
a “ fE,En * 

which comprehensively reflects the influence of the constituent stifihesw and 
the fiber volume fraction on o„. 

Consider the following composite whose fibers are aligned along the x-axis 
with E f /E m = 10, v f = 0.3 and v m = 0.35 where Ej and E m are Young’s 
moduli of the fiber and matrix, respectively, and vj and v m axe their Poisson’s 
ratios. The well-known Mori-Thnaka method (Luo and Weng, 1989) was used 


(43) 


(44) 
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to determine the five effective elastic constants of the composite strip. Figure 
9 shows a , E/E m and Eu/E m as functions of / for plane strain. Figure 10 
plots a versus / for several values of Ef f E m . Figure 11 shows the variation 
of a, E? /E m and En/E m with respect to Ej/E m for / = 0.4. It is observed 
that, when the volume fraction is low, a sharply decreases with an increase 
in /, while beyond this range a changes slowly. Besides, a decreases with 
mrrpRgiTigr Ef/E m . From the above observations it is clear that, since large 
values of A in general correspond to less stable bridged cracks, then for a 
given length of bridged crack increasing the stiffness ratio Ef/Em or volume 
fraction / will reduce the stability of the bridged crack. However, this can be 
compensated by increasing the fiber strength E, fiber radius 12 and decreasing 
the interfacial frictional friction stress r. 

4. Concluding Remarks 

A singular integral equation formulation has been presented for the anal- 
ysis of bridged cracks in orthotropic strips. The exact solutions given in 
this paper show that in terms of the normalized parameters introduced by 
Cox and Marshall (1991) the nonlinear behavior of bridged cracks can be 
depicted in a nearly universal form for materials possessing different levels 
of aniso tropy. The increase of the ratio of stiffness E//E m and volume frac- 
tion f in general reduces the stability of the bridged cracks, which can be 
remedied by increasing the fiber strength and fiber radius and decreasing the 
interfacial frictional friction stress. 
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Appendix 


In equation (9), J3>*(£) is determined by 

= /i(0> 

i=i j=i 

= /»( 0 > 

j=l i=i 

(A.1) 

+ W&«W ! ‘ = /«(€). 

jf=l i=l 

OitffcrMQ** - E^(0^ fc l = /*(0» 

i=i i=» 

where 

/i(0 = -^[Aie a ‘»« + Aae"*"*], 

/,(0 = |[A 1 Aie*»"* + A 2 A 2 e tt ’**], 

/ 3 (0 = ~~[Aie- ai( *" w)e + 

/ 4 (0 = ^[A 1 A 1 e^<'-*>* + A 2 Xae«*<*-»*]. (A.2) 
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Table Captions 


Table 1 Normalized stress intensity factor vs. orthotropy rescaling pa- 
rameters A and p for zero fiber bridged edge cracks. 

Table 2 Normalized CMOD vs. orthotropy rescaling parameters A and 
p for zero fiber bridged edge cracks. 
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Figure Captions 


Fi g. 1 Configuration, of the cracked strip. 

Fig. 2 Normalized COD for zero fiber bridged edge crack with different 
relative arthotropy. 

Fig. 3 Normalized crack month opening displacement vs. normalized 
loading parameter for bridged Griffith crack. 

Fig. 4 Normalized total stress intensity factor vs. normalized loading 
parameter for bridged Griffith crack. 

Fig. 5 Normalized shielding stress intensity factor vs. normalized load- 
ing parameter for bridged Griffith crack. 

Fig. 6 Normalized crack month opening displacement vs. normalized 
loading parameter for bridged edge crack under uniform tendon. 

Fig. 7 Normalized total stress intensity factor vs. normalized loading 
parameter for bridged edge crack under uniform tension. 

Fig. 8 Normalized shielding stress intensity factor vs. normalized load- 
ing parameter for bridged edge crack under uniform tension. 

Fig. 9 NondimensionaJ parameters a, E? fE m and En/E m vs. fiber vol- 
ume fraction f for composite with Ef/E m = 10, vj = 0.3 and 
u m = 0.35. 

Fig. 10 Nondxmensional parameter or vs. fiber volume fraction f for 
several composites with Vf = 0.3 and u m = 0.35. 
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Fig. 11 Nondimensional parameters a, E’/E m and En/E m vs. Ef/E m 
for composite with vj — 0.3, u m = 0.35 and / = 0.4. 
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Table 1 Normalized stress intensity factor vs. orihotropy rescaling parameters X 

AND p FOR ZERO FIBER BRIDGED EDGE CRACKS. 



[tt]29 , /2a 00 A 


0.355 | 0.309 
1.575 



.1 0.355 0.309 
.3 1.575 1.406 



3 

4 

0.280 

0.273 


1.286 

4.633 

4.564 

0.280 

0.273 


4.633 

4.564 

0.280 

m 


i 

4.633 

4.564 



Table 2 Normalized CMOD vs. orihotropy rescaling parameters X and p for zero 
FIBER BRIDGED EDGE CRACKS. 
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Fig. 1 Configuration of the cracked strip 




Fig. 2 Normalized COD for zero fiber bridged edge crack with different relative orthotropy 
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Fig. 9 Nondimensional parameters a, e / e,„ and E\ i / E m vs. fiber volume fraction / 

FOR COMPOSITE WITH Ef / E m = 10, vy = 0.3 AND v m= o 35 . 



f 


Fig. IONondimensional parameter a vs. fiber volume fraction / 
FOR SEVERAL COMPOSITES WITH vy =0.3 AND v m = 0.35. 





FIG. 1 1 NONDIMENSIONAL PARAMETERS a, E / E m AND E\\IE m VS. Ef / E m 
FOR COMPOSITE WITH vf = 0.3,v m = 0.35 AND / = 0.4 . 
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